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We studied a left-right symmetric model that can accommodate the neutrino dark energy {vDE) 
proposal. Type III seesaw mechanism is implemented to give masses to the neutrinos. After ex- 
plaining the model, we study the consistency of the model by minimizing the scalar potential and 
obtaining the conditions for the required vacuum expectation values of the different scalar fields. 
This model is then embedded in an SO(IO) grand unified theory and the allowed symmetry breaking 
scales are determined by the condition of the gauge coupling unification. Although 5'(7(2)_r breaking 
is required to be high, its Abelian subgroup U{1)r is broken in the TeV range, which can then give 
' the required neutrino masses and predicts new gauge bosons that could be detected at LHC. The 

I neutrino masses are studied in details in this model, which shows that at least 3 singlet fermions 

. are required. 

(N ; 

• Introduction 

c/^ : 

During the past couple of decades, astrophysical observations has improved our knowledge of cosmology tremen- 

. dously. One of the most important discovery resulting from these observations is that of the dark energy [1|. Nature 

' of the dark energy (DE) is one of the most puzzling question of physics. The observations suggest that currently , 

r~| . i.e. around redshift z ~ 1, the DE is contributing around 70% of the total energy budget of the universe, while its 

Qh" contribution was sub-dominant in the past(z » 1). Any proposed model of DE is required to satisfy these observa- 

ji^, tional constraints. These models require the mass of the scalar to be very light having scale same as Hubble scale 

^ (~ 10~^^eV). There exist myriad of such models describing the nature and the dynamics of DE (for recent reviews 

f~i , see Ref. 01] • One of the very interesting proposal for the DE is based on the fact that typical energy scale of DE 
PA '--^ (3 X 10""^ eV)'* also coincides with the neutrino mass scalepA ~ m'^. This has led to several attempts to relate 

! the origin of the dark energy with the neutrino masses [1, 0, 0, @| and this connection can have many interesting 

I> ' consequences 0,(1]. In this scenario a scalar field A called the acceleron couples with the neutrinos and consequently 

, making the neutrino mass function of A. Next, it is assumed that the dark energy poE can be written as 
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PDE = Pu + V{A). 



CNJ ' Stationary condition on pde then lead to varying the neutrino mass. These type of models are called mass varying 
neutrino (MaVaN) models 0,0, [H|- In a typical MaVaN scenario, the standard model is extended by including singlet 
right-handed neutrinos Ni, i = 1, 2, 3, and giving a Majorana mass to the neutrinos which varies with (j)a. At present 
our understanding of MaVaN models is far from being complete, several problems regarding nature origin and nature 
^ ■ of the acceleron field, about its stability etc. continue to remain. There has been a significant progress in solving 
some of these problems in the subsequent works l^, 11 1 , but much more needs to be done before this idea could be 



5_i considered as a realistic one. 
5^ , Considering the difficulties involved in constructing a reasonable MaVaNs model, most of the earlier models re- 
stricted themselves to start with the standard model and include a singlet right-handed neutrino, or else, include a 
triplet Higgs scalar. Some time back we constructed a left-right symmetric model with right-handed neutrinos and 
type-Ill seesaw neutrino masses, which could explain the dark energy with MaVaNs fl3|. In this article we work out 
some of the details of that model and embed the model in a grand unified theory. The most important feature of this 
model is that the model justifies the smallness of the very low scale, entering in this model. We have analyzed the 
consistency of the problem by minimizing the scalar potential and then have found the conditions for the required 
minima that explains the required mass scales in this MaVaNs model. We also study the gauge coupling unification 
in the SO(IO) GUT, in which this model has been embedded. The neutrino masses have also been studied and some 
conditions on the number of the singlet fermions have been worked out. 



The Model 



One of the problems with the original MaVaNs is that the condition from naturalness requires the Majorana masses 
of the right-handed neutrinos, which varies with the acceleron field, to be in the range of eV. In such models of type 
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I seesaw, the model becomes void of any seesaw, since the smallness of the neutrino masses can not be attributed to 
any large lepton number violating scale. Another restrictions of this model is that the model cannot be embedded 
in any left-right symmetric extension of the standard model, because the equal treatment of the left-handed and 
the right-handed fields would imply that if neutrino masses vary with the value of the acceleron field, the charged 
fermion masses would also vary and that would relate the scale of dark energy to the top quark mass scale, which is 
unacceptable. Although the constraint from the naturalness condition can be softened in the vDEmodels with triplet 
Higgs scalars this cannot be embedded in a left-right symmetric model. We consider here a left-right symmetric 
model, where the neutrino masses originate from double seesaw or type III seesaw mechanism and then show how 
this model can be embedded in a grand unified theory. 

In the left-right symmetric models, the standard model gauge group is extended to a left-right symmetric gauge 
group Glr = SU{3)c x SU{2)l x SU{2)r x U{1)b-l, so that the electric charge is defined in terms of the 
generators of the group as: 

Q^TsL + TsR + ^-j^^nL+Y. (1) 
The quarks and leptons transform under the left-right symmetric gauge group as: 

= ^[3,2,1,1] , Q,= (^-).[3,l,2,i], 

^.= (:^)- [1,2,1,-1] , [1,1, 2,-1], 

Sr= [1,1,1,0]. (2) 

The right-handed neutrinos Nr, is present in all the left-right symmetric model, which is dictated by the structure of the 
fermion representations and the gauge group. However, in models with type III seesaw mechanism for neutrino masses 
one introduces an additional singlet fermion Sr. As the name left-right symmetric model, the model Lagrangian is 
invariant under the left-right parity transformation given as: 

SU{2)l ^ SU{2)r, 
Ql ^ Qb. , 
h ^ in- 

In this model we have introduced the singlet field Sr, but there is no Sl. But still this model is consistent with 
left-right parity operation, since the field Sr transform to its CP conjugate state under the left-right parity as: 
Sr ^ S'^L- This also ensures that the Majorana mass term is invariant under the parity transformation, because this 
field Sr transform under the transformation SU(2)l ^ SU{2)r to itself Sr = (1, 1, 1,0) ^ (1, 1, 1,0). 
The gauge boson (excluding gluons) sector consist of two triplet and one singlet as : 



Wl, \ ( 

W^,L= I Wl^ =(1,3,1,0), W^R^ I =(1,1,3,0), S^(B-L) = (1,1,1,0) 

WZ, J \ Wr^ 

There exists several choices of the Higgs scalars, and hence, the choices of symmetry breaking chain. In the present 
model, the content of the Higgs sector will be chosen according to the following desired symmetry breaking pattern [r3|: 

SU{i), X SU{2)l X SU{2)r X U{1)^^b-l) [G3221D] 
SU{-i)c X SU(2)l X U{1)r X C/(1)(B-L) [G3211] 
"-^ SU{3)c X SU{2)l X U{1)y [G321] 

SU{3)c X U{1)q [Gem] . 

Breaking of the left-right symmetric group to G3211 requires a right triplet Higgs scalars A_r, transforming as Ar = 
(1,1,3,0). The triplet does not change the rank of the gauge group and only breaks SU{2)r U{1)r. Since it 
does not carry any U{1)b-l quantum number, it cannot give any Majorana masses to the neutral fermions. For 
the next symmetry breaking stage, U{1)r x U{1)b-l U{1)y, we introduce an SU{2)r doublet Higgs scalar field 
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Xr = (1,1,2,1) (lil . [lB |. The vev of X-R could also break [G3221_d] — > [G321], if the field Ar were not present. 
The left-right parity would then require the existence of the fields = (1,3,1,0) and xl = (1,2,1,1). Finally, 
the standard model symmetry breaking is mediated by a bi-doublet field $ = (1,2,2,0), like in any other left-right 
symmetric model. This field has the Yukawa interaction with the standard model fermions and provide Dirac masses 
to all of them. We shall introduce one more Higgs bi-doublet scalar 5" = (1, 2, 2, 0) that is needed for the purpose of 
our model. We also introduce another singlet scalar field 77 = (1, 1, 1, 0), which acquires a tiny vev of the order of the 
light neutrino masses and generate the mass scale for the dark energy naturally. 

Now we write down the explicit forms of all the scalar fields in terms of their components as 



^ - i AZ -AO ; ' ^« - V A^ -A% 



^2 n 



The most general scalar potential has to be constructed in such a way that they respect the left-right parity transfor- 
mation of the scalar fields listed below: 

Xl^Xr , Al^ Ar 
$ ^ (j)t ^ \j/ ^ \i>t 

77 <-!■ ?7 . 

Under the left-right gauge group transformation, the Higgs fields transform as 

Al ^ Ul Al Ul , Ar ^ Ur Ar Ul 

XL -^Ul XL , Xr^ Ur Xr 

In order to write down the scalar potential we also construct the fields T^$*r^ and t^"^*t^ from $ and 4" which 
transform in the same ways as $ and For convenience, we represent $ as r^$r^ as (j)2 (and similarly for \1/) 
from now on. 



Potential Minimization 



We first write down the most general renormalizable gauge invariant scalar potential respecting left-right parity 
and study details of potential minimization. Besides left-right parity, we impose following Z4 symmetry on only the 
Higgs potential to avoid few undesired terms 

XL iXL , XR -iXR , 
Al^-Al, Ar^-Ar, 



We write the the Higgs potential as a sum of of various parts and write down each part separately as: 

i,j i,j,k,l 
i,j,k,l 
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i,j,k,l 

Va = -^[tr{ALAL)+tr{ARAii)] + ^[tr{ALALf + triARARf] 
+^ [tr{ALALALAL) + tr{ARARARAR)] 

+^ MAlAl) HArAr)] 

+Y ixixL XrXr] 
Va4,4, = Yl^'l'io [AlAl + ArAr] tr{4(t>j) 

+ [^lAl + ArAr] tr{i})li)j) 

+ E f^<t>ij [ ir-(AiAi</.i,^t) + tr{ARAR4>\4>j)\ 

+ E /^■'^ii [ tr(Ai AiViV'] + tr{ARARi)\i)j)\ 

hi 

+ E'^A0y tr{AL(piAR4)'j) + E hAi;ijtr{ALilJiARipl) 

= Y.fi<l>xij [XlXl + XrXr] tr{4<Pj) 

+ E ^-/-xy IXlXl + XrXr] triiplipj) 

hi 

{Kx [xixL + XrXr] + hr,A MAlAl) + fr(AflAfl)]) rj^ 

+ E '^TO ^ [i^C-^l AlV-j) + ir(.^iAflV]) + h.c\ 
+ E [x\<PiXR + ^-c] • 



5 



We parametrize the true minima of the potential by giving vacuum expectation values to diflerent scalar fields as 
follows. 



V 0\ ^ f v' 
v' h^'^[o V 
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w J 



r/ = u . 

Since the phenomenological consistency requires v ^ w'and w ^ w' , we ignore potential terms involving w'and w' 
and write down the general scalar potential in terms of vacuum expectation values of different scalar fields 



V — V H — - V w H w 

2 4 2 4 

/^A / 2 , 2 \ , -^A / 4 , 4\ 

— Y ("i + "«) + X '"^ 

H ^- U -\ U 

2 4 

-f + + + + f (-i 4) 

+hjj u{ul + uii)vw + u{vlVr)v . 

For convenience, we have replaced + — > A^, + ^ A^, Aa + Aa Aa- The minimization of the 
potential is studied by taking partial derivatives with respect to vevs of all Higgs fields and then separately equating 
them to zero. Solving all such equations will provide us the desired values. One of the minimization conditions 

vl (^7^) ~ ''^R. (^^) ~ ^ leads to the following relation between vl and vr: 

(vr - '"D [iK ~ 9x)'"LVR - h^uv] = . 
Since (u^ — ^l) is not desirable phenomenologically, we chose 

[K~9x) 

Using above relation in an another minimization condition vl (^j^^ + ^R (^J^) = 0, we get 

4 + 4 = -^- (5) 

^x 

Parametrizing vl = A sin^, vr = A cos 9 and putting them in the two equations [4] and [5] , we find A = — /i^/A^ 
sin 26 = 29 = (f^l™) since fi^ is a large number compared to the numerator. So we get 



VR = A = y-fil/X^ , 
VL = AO ^ 



X^hx uvvr 



2 



(.9x -K) 

We have chosen the parametrization of vl and Wfl, in such a way that vr gets value equal to breaking scale of G'3211 
and Vl gets a very small value. We could have done other way around but that is not what is phenomenologically 
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allowed. Proceeding with the same kind of analysis for ul and ur, i.e., using two minimization conditions 




{S^) = and + UR (IJ^) = 0, we get 

UR = ^-Ma/^a , 

_ Aa^a (^A0W^ + hA^w^)uR 
(ffA - Aa) Ma 

Now using equation m the rj field can be shown to get vev only by term h,^u{uL + ur) as only this term is linear 
in u. The term h^u(vLVR)v does not remain linear in u after we substitute the value of vlVr from equation |4l Since 
the mass term for rj field is large and positive, we expect very small vev. So we can ignore some of the terms in the 
potential while solving for u and can easily obtain 

_ hrfVw{uL + Ur) 

M2 - (/i,,aMa/-^a) - {hr,^nl/X^) ■ 
After analyzing the complete scalar potential, we find a consistent solution with ordering 

w_R > w_R > t; > w > u > . (6) 

At this stage we can assume the different mass scales to explain the model. However, when we embed this model in an 
50(10) grand unified theory, the gauge coupling unification will impose strong constraints on the different symmetry 
breaking scales. The left-right parity and the SU{2)r breaking scale will come out to be above 10^^ GeV. So, we shall 
assume ur ~ 10^^ GeV. We also assume to^ ~ ttt-a ~ ur. However, it will be possible to keep the G3211 symmetry 
breaking scale to be very low, and hence, we shall assume r-^ vr ^ TeV. We find the remaining mass scales to be 
V ^ rriw ~ 100 GeV, m ~ ul ~ eV and vl ~ 10~^ eV. 



Embedding The Model In 5'O(10)GUT 

The idea of Grand Unified Theories (GUTs) has emerged as a very attractive idea to go beyond Standard Model 
(SM) for last three decades. It unifies the three different looking gauge coupling constants of the SM, and in addition, 
reduces the number of particle irreducible multiplets into lesser number of multiplets. The ad-hoc looking hypercharge 
assignment in SM gets a predictive framework in GUTs, i.e, the charge quantization remains no more a surprise in 
GUTs. The smallest GUT SU{5), in its non-supersymmetric version, does not unify the three gauge coupling constants. 
Out of the higher rank gauge groups containing SM gauge group as a subgroup, the rank four semi-simple group 5*0(10) 
has emerged as a very attractive candidate for GUT. It can accommodate the entire SM fermion content in its single 
16-dimensional complex irreducible spinor representation including the right handed neutrino with three copies for 
the three families. Its all irreducible representations are anomaly free providing a natural predictive framework to 
understand the fermion masses and mixing. Also the seesaw structure gets a natural embedding in 50(10). The 
left-right symmetry group can also be embedded in S'O(IO) GUT. 

We shall study here the embedding of the present model with all its Higgs content in S'O(IO) GUT. We consider 
the following breaking pattern of S'O(IO) gauge group to first Pati-Salam gauge group 5C/(4) x SU{2)l x SU{2)r, 
next to the left-right gauge group and then to the SM gauge group 



Mu 


SU{A) X SU{2)l X SU{1)r [G422D] 




Mc 


5t/(3)c X SU{2)l X SU{2)r X C/(1)(b_l) 


[G3221D 


Mn 


5(7(3)c X SU{2)l X U{1)r x U{1)(b-l) 


[G3211] 


771 J- 


SU{i)c X SU{2)l X U{1)y 


[G321] 


raw 


5(7(3)c X U{1)q 





The Higgs multiplets which can provide the masses for all the SM fermions are limited as 16 x 16 = 10s + 120a + 126^. 
The 10 dimensional Higgs field decomposes under left-right gauge group as 

(10) = $(1, 2, 2; 0) © (3, 1, 1; -i) ® (3, 1, 1; i) . 
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One can oasily identify the bi-doublet $(1, 2, 2; 0) appearing in the left-right model contained in i?$(10). To include 
another bi-doublet \1/(1,2,2;0) present in the model, a second Higgs field i?.t(10). 

Although the fermion and gauge sector of the 50(10) GUT model are quite simple, the Higgs sector is quite 
complicated since it is not only required for generating fermion Masses, but an appropriate Higgs content is also 
needed for systematic and consistent breaking of the 50(10) gauge group down to the SM gauge group in one or more 
steps. To break 50(10) gauge group to the Pati-Salam gauge group, one requires Higgs field either 5(54) or T(210), 
which decompose under Pati-Salam group as 

5(54) = (1,1, 1)8 (1,3, 3) 8(20,1, 1)0(6,2, 2), 
T(210) = (1,1,1)8(15,1,1)0(6,2,2)0(15,3,1) 
8(15,1,3)8(10,2,2)8(10,2,2). 

Giving vev to either of the two fields in the singlet direction will serve the purpose of the desired breaking. The 
(15, 1, 1) of T also has a singlet under the left-right gauge group which can acquire vev to break the Pati-Salam group 
to the left-right group. The (15,3,1) and (15,1,3) Higgs multiplets of T also contain the fields At(l,3, 1,0) and 
A/j(l, 1, 3, 0) present the left-right model. Hoever, the T singlet under Pati-salam gauge group is odd under D-Parity. 
If we give vev to T singlet, the left-right symmetry will be broken at unification scale itself. Since our model is 
left-right symmetric, we must avoid D-parity breaking until left-right group is broken. 

However, the singlet in 5(54) field under Pati-Salam gauge group does respect and so can be used to break the GUT 
group to the Pati-Salam gauge group. But, the breaking with 5(54) does not serve the purpose of further breaking 
to the left-right group. So for the next step breaking, a Higgs Field ^4(45) is needed along with 5(54) which has the 
decomposition under the left-right group as 

A(45) = (l,l,l;0)8Ai(l,3,l;0)8AK(l,l,3;0) 
8(3,1,1;^)8(3,1,1;-^)8(8,1,1;0) 

8(3, 2, 2;^) 8 (3, 2, 2;-^). 

The first row of the above decomposition is of our interest as it contains the fields Ai:,(l, 3, 1, 0) and Ai^(l, 1, 3, 0) of 
our model along with the left-right group singlet. This singlet is even under D-parity and so the left-right symmetry is 
unbroken until A^ acquires vev along the singlet direction to the SM gauge group. We will be following this approach 
in the remaining part of this section. 

Now the fields Xl (1, 2, 2, 1) and Xni^i 1) 2, 1) are still left to be embedded in some tensors of 50(10). The desired 
quantum numbers indicate that they can be embedded in the spinorial Higgs representation (C(16) 8 0(16)) . 
Decomposition of the 16 8 16 spinor representation under left-right group are given as 

0(16) = XL(l>2,l,-l)8Xfl(l,l,2,l) 
8(3,2,1,^)8(3,1,2,-1), 

0(16) = XL(l,2,l,l)8Xfl(l,l,2,-l) 
8(3,1,2,1)8(3,2,1,-1). 

Having embedded all the Higgs fields of our model into 50(10) tensor fields, we now write vacuum expectation 
values along the three singlet direction under the SM group of the fields A(45) and 5(54) as 

(A) = McAc + MrAr, 
(5) = MuS, 

where Ac, Ar and 5 are the singlet directions under the SM gauge group given as 
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Ac = i^Ase + Ajs + Agioj 

Ar = (il2+i34) 

4 10 
= 3 X ^ ^ .Saa — 2 X ^ ^ iSaa . 
a— 1 a— 5 

The indices (1, 2, 3, 4 ) belong to S'0(4) and (5, 6, 7, 8, 9, 10) belong to S0{6) subgroup of the group 50(10). We have 
not taken care of the normalization factors while writing the directions of the singlets as they are not much relevant 
for the present discussion. However, we can assume that the normalization factors are absorbed in the corresponding 
vev values and can proceed without worrying about them for an approximate analysis. 

Let us denote — h, i/* = H for simpHcity in notations. Now we write the most general S'O(IO) invariant Higgs 
potential: 



V = pi\ AabAba + Ms SabSba + ^^l haha + fJ-H HaHa + (CC) + /l^ 7?^ + 'q^ 

+ Xa A^A^ + A^A^ + \sS^ + \h h* + XhH^ + A, {CCf + A^ {C* + C*) 

+ OAsA^S^ + g'AsAabAbcScdSda + QAsAabSbcAcdSda 

+ ha {ghAAabAbc + ghsSabSbc) he + {g'hAA^ + g'hsS^) h^ 

+ Ha [gHAAabAbc + gHsSabSbc) He + [gH aA^ + 5hs5'^) H"^ 

+ {ghch^ + gncH^ + gAcA^ + gscS'') CC + g^HCn h {CC + CC) . 

The ^4 symmetry (expression[3l) used while writing the Higgs potential invariant under left-right gauge group has also 
been imposed here on the corresponding S'O(IO) Higgs multiplets. Moreover, we have prevented some of the terms 
by applying the discrete symmetry S ~^ ~S. The realization of the first three symmetry breaking steps is possible by 
taking the following structure of the vev assignments to the fields ^(45) and S'(54): 



(A) = iT2 (E) diag{MR, Mr, Mc, Mc, Mc) 

/ 3 3 
(S) = I®diag{--Mu, --Mu, Mu, Mu, Mu 

For the matter of convenience we have just replaced the weus with the corresponding breaking scales. The potential , 
in terms of the vev values of A and S, will be approximately given as 



V = ii\ (6A4 + 4Af2 ) + ^2 i5^^2 ^ ^^2^ ^ gAcQMl + gsc ISM^ ) CC 

+ (m^ + ghS^Ml) haha (a = 1 - 4) + {nl + guAQMl + ghs^M^) haha (a = 5 - 10) 
+ [A + gns^Ml) HaHa {a^l-A) + {^?H+ gnAGMc + gHS&Ml) HaHa (a = 5 - 10) 
+ A^ {Ull + AMlf + X'a (6M4 + 4M4) + XgM^ + gAsMfj (6Mg + 9Af^) 
+ Xhh^ + XhH^ +gr,HcTlh {CC + CC)+Xc {CCf + X'^ (C^ + C^) + A^ 77^ 

We have assumed Mr_ <^ Mu ^ Mc while writing the final form of the potential. In order to give desired masses 
(of the order of Mw) to the two left-right bi-doublets , fih and /i// will have to be fine-tuned at the order of 
scale of Mjj- The fine-tuning can produce very large masses to the triplets of /i( or H) provided the condition 
(^gfiA^M^ — hhsiM^^ ~ (-|-M^) is satisfied. Another fine-tuning is required in the mass parameter Mc to provide 
the desired TeV scale masses to the Higgs fields C(16) ® C(16). Before ending this section, we would like to notice 
an important point. If we take the gsA coupling to be very small, we can argue that the appearance of the similar 
combination (6Mp -I- 4M^) everywhere in the potential allows Mc and Mr to take quite different values without 
disturbing other part of the potential. So the scale of Mc and Mr can be chosen to be different by orders of magnitude 
to get the desirable breaking. 
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Gauge Coupling Evolution 



In the present section, we will be studying the set of two- loop renormalization group (RG) equations for the evolution 
of the coupHng constants and will be verifying the consistency of the chosen vev for different Higgs fields in the context 
of 5*0(10) GUT. For simpHcity, we assume that the scale Mu and Mc are very close and we ignore the evolution of 
the coupHng constants between the two scales. This is quite preferable as we will see later that the unification scale 
is very tightly constrained by the current proton decay bound and any substantial difference between the two 
breaking scales would make it even worse. We start with the following equation for the two-loop evaluation of the 
coupling constant Ui 




dt I _-l I 

where t — In (Af^) and is the desired energy scale where the coupHngs constants, a^'s, are be determined. The a^'s 
and bij's are the one-loop and two-loop beta functions governing the evolution of a^'s and include the contributions 
from gauge bosons, fermions and scalars in the model. 

The fermion contribution to the beta function is taken right from the starting, the electroweak scale (lOOGeV). The 
contributions of the gauge bosons to beta functions are straightforward to compute as one can easily determine the 
expected mass scales of the heavy gauge bosons corresponding to any given gauge group. However, the contribution 
coming from the Higgs content is not so clear because the heavy Higgs modes can have various possible mass spectrums. 
We will use the extended survival hypothesis to fix this uncertainty. The extended survival hypothesis is based on 
the assumption that only minimal number of fine-tunings of the parameters in the Higgs potential are imposed to 
ensure the hierarchy in various gauge boson masses. According to the extended survival hypothesis, only those scalar 
multiplets are present at any given intermediate breaking scale M/ of a intermediate gauge group Gj which are either 
required for breaking the gauge group Gj or needed to further break any other intermediate gauge group below scale 
Mi. 

A list of Higgs multiplets surviving at the breaking scale of a intermediate group G/, using the extended survival hy- 
pothesis, are given in table. A list of both one-loop and two-loop beta coefficients, which include all the contributions, 
that govern the evolution above the breaking scale of G/ to the next intermediate scale are also listed. 

Since our model contains intermediate steps, we require appropriate matching conditions at the corresponding 
breaking scales. For the tow- loop RG running of the coupling constants, the matching conditions have been derived 



m 



isl . lil . Suppose a gauge group G is spontaneously broken into a sub-group Yii Gi with several individual factors 



Gi, then the following matching condition need to be satisfied for the two-loop analysis 



where C{G/Gi) is the quadratic Casimir invariant for the group G/Gi. We choose initial starting values of the 
above three coupling constants ( central values) at scale Mw to be a^y {^w) — 59.38, a^2(-^w) — 29.93, and 
a'^^iyMw) = 8.47. Now let us write the 

The boundary conditions at various breaking scales, using the expression [HI can be written as 

1. At scale m^: 

3 2 

2. At scale Mr: 



^ik (Mr) = «2K (Mr) ~ , 
"2"fl (Mr) = {Mr) . 
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Group Gi 



Higgs content 



1) io ® (1' 2' -|)io 

2) 10' ® 2' ^2)10' 
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V ' J 



/ -26 I ii 
12 



8 I 



\ 44 18 104 / 
^ 5 5 25 ^ 



G3 



(1,2, i0)^„©(l,2, , 
(1, 2, ,10) ©(1,2 

[1. 1. -il)i6 + (l^l ' 



' 2 ' 2 J 16 



-3 



53 
12 



/-26 I 
12 8 



12 

4 



o ir 15 

A Q 



15 65 / 



G3221 



(1, 2, 2, 0)^0 
(1, 2, 2, 0)10, 
(1,2, . 
(1> 1. 2, i,^^ 
(1, 1, 3, 0)45 
(1, 3, 1, 0)45 



1' -l)i6' 



HI' 2, L 
(1, 1, 2, 



' 27T6 



-26 
12 
12 
4 



9 

2 


9 
2 


1 
2 


\ 


39 


3 


9 




2 




4 




3 


39 


9 






2 


4 




27 


27 


23 




4 


4 


4 


/ 









Table I: Higgs multiplets at different intermediate breaking scales along with the both one-loop and two-loop beta coefRcientss, 
including all the contributions from fermions, gauge bosons and Higgs bosons, which govern the evolution of coupling constants 
above breaking scale of G/ to the next breakingscale. 



3. At the unification scale Mjj 
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The matching conditions at the unification scale have been written by assuming the Pati-Salam scale to be almost 
close to the unification scale. 

Using the above boundary conditions we have numerically solved the equation[7]for the two-loop RG evolution for all 
the coupling constants. We have taken the breaking scale of the gauge group G3211 to be around ITeV. The unification 
scale comes out to be Mu = lO^^ ^'GeV and the corresponding coupling constant is estimated as {Mu) = 43.4. 
Also the breaking scale of left-right symmetric gauge group, i.e., 6*32211) turns out to be M^ = 10^^'^GeV. The running 
of the various coupling constants with energy scale are shown in figure [H 

However, we find that the scale of the unification along with the — 1 are not satisfying the most recent bounds 
on proton decay, although very close to the limit. The current experimental lower bound of the partial life time for 
p e+7r° is Tp > 8.2 x 10'^^ years and for p /i+Tr" is Tp > 6.6 x 10^'^ years The theoretical decay rate of the 
proton can be estimated as: 
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Figure 1: Evolution of coupling constants 
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Table II: Threshold contribution at left-right breaking scale 



This can be used to estimate the lower hmit of the Heavy gauge boson masses. If the mass scale of super heavy 
^auge bosons are given as Mx — 10"GeV, the above proton decay bound is equivalent to 



/ aGUT \ 

V 45 / 



X 



^q2(«-15) > 8^ (9) 



What we obtain for the value of n in our analysis is k = 6.07. This is below the lower limit allowed by the proton 
decay bound as specified in the right-hand side of the expression [9l However, the value of k is very close to the 
allowed lower limit and so we will try to explore the viability of our model by allowing threshold uncertainty in the 
Higgs spectrum at various intermediate breaking scales. It is important to remark at this point that we could get the 
reported value of k to be close to the limit only when we optimized certain degrees of freedom in the Higgs sector. 
For instance, the Higgs-bidoublet <& has been asumed to arise from a real 10-dimensional 5*0(10) Higgs H^. So <& 
would not be equivalent to two SM Higgs doublets at the electroweak scale but will be equivalent to only one such 
doublet. Similar asuumption has been also taken for \I>. However, we would Hke to emphasize that the results and 
discussion of the potential minimization will remain almost same. 

The threshold uncertainty in the Higgs spectrum arises form the fact that the Higgs bosons becoming heavy at 
a given breaking scale may not get exactly same masses equal to the energy corresponding to the breaking scale. 
However, the Higgs mass spectrum is expected to be scattered around the energy of the breaking scale within an 



small width. For our analysis, we follow a similar approach discussed in [2l|. We assume that the masses of the Higgs 
bosons are scattered around the breaking scale within the factor of ^ to 30. So if the mass of a Higgs multiplet at 
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Table III: Threshold contribution at the unification scale 



the given breaking scale Mi is Mh, then we expect 



1<^<30. 
30 Mi ^ 

To include the threshold uncertainty at a given breaking scale, we need to slightly modify our matching conditions 
at that scale. The matching condition given in expression [8] is modified as 

where Xi = ajn^^^. So the threshold uncertainty has been included in the matching condition due to presence of the 
term involving ln(M/f/M/). 

To avoid any over estimation of the threshold uncertainty we assume that all the Higgs multiplets, belonging to a 
single common irreducible Higgs representation of S'O(IO), becoming heavy at a given breaking scale will have the 
same mass scale around the breaking scale. 

The threshold uncertainty at the breaking scale of gauge group G'32iiis vanishing. The Higgs multiplets, coming 
from different 5*0(10) irreducible Higgs, contributing to the threshold uncertainty at remaining two intermediate 
scales, the left-right breaking scale and the unification scale, are listed in the table [ll] and IHIj respectively. The 
corresponding calculated beta-coefficents, (ai)'s, which include the contribution from all the Higgs multiplets coming 
from the same 50(10) irreducible representation (as their masses are assumed to be same), are also shown for the 
two breaking scales. 

Now using these calculated a^'s and including uncertainty in Mh/Mi, as discussed before, we have shown a scatterd- 
plot between coupling constant a^^ and the corresponding unification scale Mu in figure [2j We have numerically 
obtained the values for and Mu for randomally chosen values for Mh/Mi between the range (^ — 30). The 
random values for all the Higgs multiplets belonging to the same S'O(IG) ireducible Higgs are taken to be same at 
one perticular breaking scale but different at the other breaking scale. 

Moreover, we have aslo plotted the curve corersponding to the most recent proton decay bound (red solid curve) 
[l3| and relatively older proton decay bound (blue dashed curve) [i^l in figure[2]to show the allowed region in a^^-Mjj 
plane. Only the right part of the curve is allowed by the bound. It is worth noting that the allowed parameter space 
is more and more constrained as more updated data on proton decay bound is available. However, we get a resonable 
allowed region in the figure [21 although small, even after allowing the most conservative threshold uncertainty. So 
we expect our model to be satisfactory within the tolerable amount of threshold uncertainty as far as proton decay 
bound is concerned. 
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Figure 2: Threshold uncertainty in the unificaton scale. 



Yukawa Sector And Neutrino Masses 



In the present section, we discuss the origin of neutrino masses in the model. Before proceding further we would like 
to make it clear that the discussion about neutrino masses in the present section will only move around the left-right 
symmetric model with few inputs from the S'O(IO) GUT in motivating about certain patterns for taken Dirac mass 
matrices for fermions in our analysis. Moreover, the discussion will be mainly focused on the matrix structure of low 
energy neurino mass matrix allowed with certain assumptions. We will aslo argue, in what follows, that the consistent 
neutrino mass spectrum is not possible within picture of one or two SO(IO) singlet fermions S. We start by writing 
the Yukawa sector of the model as 

Cy = r,, Ir,^ + r/, JI:^ iRj^ + {Fl),„ IuXl + {Fr),^ S^Jm XR (10) 

+ ^M,nn1lS''LmSRn (11) 

The Yukawa coupHngs Y and Y' are 3x3 matrix, while Fl and Fr are 3 x n matrices, if we assume that there are 
n singlet fermions S. So M is a n x n matrix. Our study of consistent embedding of the model in SO{10) GUT 
requires same structure for both Fl and Fn up to the scale of left-right symmetry breaking which, after RG running, 
can produce small difference at the weak scale. For the present discussion we assume it to be small enough so that it 
can be safely ignored. 

The Dirac masses for all the SM fermions including neutrinos are generated form the the first two terms by giving 
vev to the bi-doublets as in any other left-right symmetric model. Since $ and 5* are coming from two independent 
and real 50(10) 10-dimensional Higgs, the Dirac mass matrix for neutrinos and charged leptons are independent. 
However, the Dirac mass matrix for the up-type quarks have the same structure as the Dirac mass matrix for the 
neutrinos and similarily the Dirac mass matrix for the down-type quarks will have similar structure as the Dirac 
mass matrix for the charged leptons (simply because all SM fermions are assigned to a multiplet of S'O(IO) GUT). 
Although, these similarities in the structures are exact only at the GUT scale, we expect some of its features to be 
more or less same even at the low scale. So we can well assume that the Dirac mass matrix of the neutrinos would 
almost appear diagonal in the basis where the charged lepton mass matrix is diagonal. The assumption is based 
on the observation that the up-type and down-type quarks are simultaneously diagonal in the a basis as the quark 
mixing matrix is very close to unity. So we borrow the pattern from the quark sector to the lepton sector where the 
structure of Dirac mass matrix of the neutrinos is not directly known unless neutrinos are Dirac fermions. We expect 
the following pattern of the Dirac mass matrix of neutrinos in the diagonal basis of the charged leptons 
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M^D = vYi 



lepton 



nit 



rrib 
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mr 




^ 0.0001 
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y 0.3 



where and mi, are masses of top and bottom quarks and mg, m^, to,- are masses of electron, muon and tau leptons. 
The part of the Lagrangian relevant for the neutrino mass generation is given as follows, 



s)^ 



iVf, S„ 



X 




H.C. 
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V -^mj^i Pmj'^R M^nU 







FinVR 



TV? 
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H.C 



(12) 



(13) 



Our first task is to analyze the mass spectrum provided by the matrix X in case of one generation of all fermions. 
We write the eigenvalue equation as (eigenvalue: A): 



- Mu A2 - i^^wlA - 2YF^vvlVr - MY^uv^ = 

Case 1: A >> v, we get 

X{X + FvB.) {X- Fvr) = 

The above eigenvalue equation predicts two TeV scale Majorana fermions. The massless solution contradicts with the 
condition we started with, and so is unphysical. 



Case 2: A << v, we get 



2YvvL MY^uv^ , 
^ = + Z.2 2 14 



which is of order of eV. So the two Majorana fermions pick up masses of the order as high as TeV and one remains 
sufficiently Hght (~ eV) to be identified as light neutrino. 

To make the discussion some more general, we take three generations for all the SM fermions including the left and 
right handed neutrinos but only one generation for the singlet S. We look for a possibility whether it can account 
for the existing picture of three light active neutrinos. To search for any such possibility, we try to find out the mass 
spectrum, within this scenario, by solving for the eigenvalues of the matrix X. To simplify further, we take all the 
eigenvalues of the matrix M^d to be same with a common value equal to the largest one for initial analysis. This 
enable us to factor out (A^ — z'^v'^) from the algebraic expression of Det {X) predicting four Majorana fermions of 
scale around 10 GeV. The rest of the factors have got the same form as the expression of determinant in case of one 
generation of all SM fermions, as discussed earlier, leading to the two TeV and one eV scale Majorana fermions. The 
scenario provides us only one light neutrino and, hence, can not account for the observed neutrino mass spectrum. To 
explore the effect of some possible hierarchy present in the eigenvalues of the Dirac mass matrix of the neutrino like 
one present in the charged lepton mass matrix, we take two of the eigenvalues to be same and vary their scale below 
the third one. We are still able to expHcitly get two of the Majorana fermions having mass scale equal to me j • 
One may think that the remaining two Majorana fermions might get mass scale as light as eV leading to three light 
neutrinos. To rule out any such possibility, we have plotted the masses of the two remaining Majorana fermions 
(which comes out to be same) with the ratio of the two mass scales of the eigenvalues of the Dirac mass mass matrix 
of the neutrinos in figure [31 We find that the masses do not go below the lightest mass scale of the eigenvalues of 
TOi/D. Even in two generation scenario of S fermions, there is not much progress except we get two eV scale Majorana 
fermions which is still not sufficient. 

We now turn to the case of three generation for S fermions. One obviously expects to get the three light neutrinos. 
The basic way to get the low energy neutrino mass matrix has been outlined in [l3| which is given as 
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(r + Y'^) + rY {FM~^F^y^ 



(15) 



as we have uv'^ — r vvlvr in our model (expression 4]) where r = (A^ — g^) /h^. 

The first term is the type-Ill seesaw contribution 23| and the second term is the double seesaw contribution. With 
the choice of the vevs, it is obvious that this scenario provides us with three eV neutrinos. 

Now we will try to explore the limits of the expressionfTSlfor low energy neutrino mass matrix to check its consitency 
with current data on neutrino masses and mixing by allowing some very simple form for matrix M. In the basis where 
charged lepton mass matrix is diagonal, the neutrino mixing matrix {Upmns) is just the matrix that diagonaHzes the 



{UpMNsf^ rriuUpMNS = mf"^ = Diag {mi,m2,mz) 
The UpMNS mixing matrix is usually parametrized in the literature as 



UpMNS = i?23 (^23) Ri3 iOi3,S) i?i2 (^12) -Dag (e^"! , e*"^ , l) , 

where Rij are the rotation matrices in the ij plane with angle 9ij. S is the CP violating phase associated with 1-3 
rotation and 77's are the Majorana phases appearing only in the case of Majorana neutrinos. To date, two mass square 
differences and three angle have been measured while CP violation is completely unknown in the leptonic sector. We 
take the following observed values for three mixing angle and two mass square differences at 90% confidence level 
from particle data group [2J| as: 



Am^i^m^-mf = (8.0 ± 0.3) x 10"^ eV^ 
Aml^ ^ml-mj = 1.9 to 3.0 x 10"^ eV^ 
sin2(20i2) = OMtofA 

sin^ {023) > 0.92 

^(6*13) < 0.19 

We will be mainly using the mean values of the observed parameters in our analysis. 
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In its most general form, it is straight forward to argue that m„ can accommodate the existing data on neutrino 
masses and mixing simply due to the presence of enough number parameters in F and M unless type III term 
dominates significantly. An interesting thing would be to consider some simpler form of the neutrino mass matrix 
by reducing appropriate number of parameters with some tolerable assumptions. The basic idea is to explore the 
possibility of any such simpler structure in light of the current neutrino oscillation data. 

We start with the assumption that the three singlet fermions S are blind to their generation within themselves 
leading to the following democratic structure of matrix M : 




M = 

The structure allows us to believe that there is no induced mixing between the left-right neutrinos and the singlets. 
So, F matrix can be written as product of a unitary matrix and a diagonal matrix. The unitary matrix connects the 
basis of the democratic structure to the basis where the charged lepton mass matrix becomes diagonal. To get some 
more simpHcity, we are driven to assume that the two basis are identical, i.e., the unitary mass matrix is identity 
matrix. It leads to the following structure of the low energy neutrino mass matrix: 



nit , 
nib ' 



where a, (3 and 7 are the final parameters appearing in the neutrino mass matrix after absorbing all the parameters 
present in F, M and Y . We take the following familiar tri-bimaximal form of |2H| of the Upmns mixing matrix for 
our discussion and attempt to diagonaHze m,^ having above structure: 



UpMNS — Utbm — I —1 



V2 









-V2 


V3 



where 6*23 — 7r/4, 6*13 — 0, and sin^ 612 — 1/3. 

We attempt to diagonalize m^, with the tri-biamaximal form of the mixing matrix which requires the following 
relation of the parameters a, /? and 7 with masses of the charged leptons as: 



a 



(3 ~ ^ ~ 0.05 



-0.75 



The diagonal neutrino mass matrix comes out of the form: 



Daig 



2mt 
nib 



( rUe 



V 





me 
2 








/ vvl_ 



So the present form of and U pmns produces degenerate masses for the two light neutrinos which is likely to be 
cured once we sHghtly deviate from tri-bimaximal form of Upmns- The deviation can be realized either by taking 
non-maximal value of 6*23 or non vanishing value of 6*13 or both. We take only non-zero value of 6*13 to be the sole 
realization of the deviation for our purpose. The deviated form of tri-baimaximal matrix for very small value of ^13 
can be parametrized as: 
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^ V l-%/20i3 -V2-013 V3/ 

While trying to diagonalize the m^, numerical methods are used to find out the desired values of the free parameters. 
We find that the degeneracy encountered in the case of tri-bimaximal mixing matrix disappears as soon as finite value 
of is introduced. This finite value is determined by imposing the condition Am^i/Am^i ~ 0.033 which leads to 
following value of sin 613 

sin 6*13 = 0.11. 

The value is well within the allowed value for ^13 from oscillation data. The correct scale of the mass square differences 
is easily achieved by adjusting the over all scale of the neutrino mass matrix. The corresponding values of the other 
parameters come out to be 

a = 0.02 
(3 = 0.06 
7 = -0.75 

The point we would like to emphasize is that even the simple structure of the mass matrix taken in our analysis 
is able to account for the existing framework of three active light neutrinos even though the assumptions may not 
correspond to any real underlying symmetry. 

Dark Energy 

We shall now show how the model can accommodate the proposal of the mass varying neutrinos (MaVaNs) 0, Q. 
The basic idea behind the mass varying neutrinos is that some scalar field, the acceleron, acquires a value of the order 
of 10"'^ eV, which gives an effective potential that contributes to the dark energy with the equation of state uj = —1. 
However, till recently the neutrino masses were contributing to the effective potential much more strongly and the 
combined fiuid of the background neutrinos and the accelerons were behaving as dark matter with the equation of 
state Lu — 0. As the neutrino masses were varying with time, the contribution of the background neutrino density to 
the effective potential were changing. Only in the recent past, the contribution of the acceleron field to the effective 
potential became stronger than the background neutrinos, changing the equation of state of the combined fluid, and 
the universe started accelerating with dark energy domination. This can then explain why the scale associated with 
the amount of dark energy is comparable to the neutrino masses, why the amount of dark energy is comparable to 
the ordinary matter, and why the universe is dominated by dark energy only now and for the rest of the time in the 
past the evolution of the universe was governed by matter. 

In spite of these advantages, the MaVaNs scenario are not free of problems. We shall now try to explain how 
the MaVaNs scenario can be accommodated in a grand unifled theory. After describing the generic features of the 
MaVaNs, following the original proposal [4], we shall explain how our present model answer this question. We shall 
not restrict ourselves to any particular choice for the acceleron field, and hence, consider the potential for the acceleron 
field to be same as that considered in the original proposal. At the end we shall mention how the present model can 
be extended to allow a milli-eV mass pseudo-Nambu-Goldstone Boson (pNGB), which can become the acceleron field. 

We shall now mention this possible origin of the acceleron field in an extension of our model. Following the 
prescription followed in ref. 0, we introduce three 77 and several Higgs doublets. The vevs of the fields 77 would then 
give rise to global symmetries, which are allowed by all the Yukawa couplings due to the choice of quantum numbers 
of the Higgs doublets under these global symmetries. However, when the Higgs doublets acquire vevs, the global 
symmetries will be broken and there will be pNGBs, which couple to the neutrino masses. Although the dynamics of 
the pNGBs are not specified, the masses and the potentials of the pNGBs are determined by the Coleman- Weinberg 
potential, as demonstrated in ref. [^. Since the introduction of several Higgs doublets may not allow the the gauge 
coupling unification, we shall not discuss this extension any further. Moreover, there could be some other origin of 
the acceleron field, so from now on we shall only mention the generic features of this model. 
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In a generic MaVaNs models, the coupling between neutrino mass and A induces the following effective potential 

V = (p, - 3P,) + Vo{m,) (16) 

Here the scalar potential Vo{m,j) is due to the acceleron field (written as a function of neutrino mass) and Pi, is 
pressure of the neutrino fluid. In the late time evolution the non-relativistic Hmit i.e. nii, T is of particular 
interest. In this case Pu ^ and one can write the effective potential as, 

V = mi,n„ + Va{m,y) . (17) 

The acceleron field will be trapped at the minima of the potential, which ensures that as the neutrino mass varies, the 
value of the acceleron field will track the varying neutrino mass. One can write equation of state in the non-relativistic 
case for a combined fiuid of neutrino + acceleron; 

w = P/p = — (18) 

rrii^rii, + p^ 

One generic feature of this solution is that it gives w w —1 at present. The most important feature of this scenario 
is that the energy scale for the dark energy gets related to the neutrino mass, which is highly desirable. This also 
explains why the universe enters an accelerating phase now (j^ . 

We shall now discuss the implementation of the z^DE mechanism in our model. For simplicity, we consider only 
one-generation scenario. The effective scalar field potential of the scalar is of the Coleman-Weinberg type i.e. 

Vo^AUog{l + \MM)/f^\ (19) 
where, Mg is the singlet fermion mass. We assume that Ms{A)/p. 3> 1. 

Ms = M (77) = M u 

depends on the acceleron field A. Thus the neutrino mass becomes a dynamical quantity. When the neutrinos become 
non-relativistic the dependence of Mg on A governs the dynamics of the dark energy. A is chosen in such a way to 
ield the dark energy density JIde ~ 0.7. This type of potentials are extensively used in the dark energy literature 



s 



27] . Now we can write the effective low-energy Lagrangian in our model 



y2 ^2 

Cejf = Ah{A)-^— v,v^ + B.C. + A4log(l + \Ms{A)IIl\) , (20) 



From the choices we have made about the vevs, we have retained only the dominant double seesaw term [14] in 
the effective Lagrangian. As u 0{eV), the mass parameter is of the order of eV. Since the ratio {v/vr)^ ~ 
10~^ — 10"'^, the Yukawa couplings coupling to be of order unity. Thus the first two terms in equation (14) are 
comparable to the last term describing the dark energy potential. 

The Majorana mass of neutrino varies with the acceleron field through the parameter Ms and the mass scale of this 
parameter remains near the scale of dark energy naturally. The interesting feature of our model is that we do not 
need any unnaturally small Yukawa couplings or symmetry breaking scale to achieve this naturalness requirement. 
Also the variation of Mg does not affect charged fermion masses in the model. Moreover, the electroweak symmetry 
breaking scale v and the J7(1)_r breaking scales are comparable and hence the new gauge boson corresponding to the 
group U{1)r will have usual mixing with Z and should be accessible at LHC. 

Since the local minimum of the potential relates the neutrino mass to a derivative of the acceleron potential, the 
value of the acceleron field gets related to the neutrino mass. The acceleron field provide an effective attractive force 
between the neutrinos. When this effective force is stronger than the gravity, perturbations in the neutrino-acceleron 
fiuid become unstable. The source of the free-energy comes from the attractive interaction between the neutrino and 
the acceleron field. The instability is similar to that of the Jeans instability found in a self-gravitating system. The 
instability can lead to inhomogeneity and structure formation; the instability would grow till the degeneracy pressure 
of the neutrinos would arrest the growth. The final state of the instability would produce neutrino lumps or nuggets 
|27, 2^. The neutrino lumps would then behave as dark matter and will not affect the dynamics of the acceleron field 
[28l |. This instability is a generic feature of MaVaNs scenario, however it can be suppressed if the neutrino become 
superfiuid [i^ or if the MaVaNs perturbations become non-adiabatic. 
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Conclusions 



We have constructed a left-right symmetric model of i/DE that can be embedded in an 50(10) GUT. After 
discussing the Higgs content needed for the model, details of potential minimization have been carried out considering 
all possible allowed terms. In particular, we have tried to explore the possibility of choosing the minima such that 
only neutral Higgs components get vev without constraining the couphngs constants. But it turns out that some such 
constraints are needed in most general form of the potential. The complete analysis allows the desired ordering of the 
vevs. Then we study the embedding of this left-right symmetric model in 5*0 (10) GUT. We show that SO (10) GUT 
with Higgs multiplets S'(54), A(45), two H{10), C(16) © C(16), 77(1) along with an additional fermion singlet is able 
to accommodate the left-right symmetric model. The embedding allows the Pati-Salam and the left-right symmetry 
group breaking scales to be different by orders of magnitudes. We have studied the one loop RG running of various 
couplings constant and have found that the desired assignment for vev values for different Higgs fields is consistent 
with the gauge unification. Then the origin and possible structure of neutrino masses and matrix have been discussed 
in detail. It has been shown that generation of three light active neutrinos of eV scale is not possible in scenario 
with one or two 5*0(10) singlets fermions. In the generic case of three singlets, we have taken a simple structure of 
neutrino mass matrix with some tolerable assumptions and shown that the structure is consistent with current data 
on neutrino masses and mixing. Then we described implementation of uDE in the model. The model allows the 
mass parameter of the singlet, which varies with the acceleron field, to have the same scale as the scale of dark energy 
satisfying the desired naturalness requirement. 
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